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Abstract. In this paper we study ensembles of random symmetric matrices 
X n = {Xij}™ J=1 with dependent entries such that KXij = 0, E_Xy = cry, 
where Oij may be different numbers. Assuming that the average of the 
normalized sums of variances in each row converges to one and Lindeberg 
condition holds we prove that the empirical spectral distribution of eigen- 
values converges to Wigner's semicircle law. 
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1. Introduction 

Let Xjk, 1 < j < k < oo, be triangular array of random variables with E Xjk = 
and El? fe = cr? fe , and let Xjk = X^j for 1 < j < k < oo. We consider the 
random matrix 

X n = {Xjk}™^ = i- 

Denote by Ai < ... < X n eigenvalues of matrix n _1 / 2 X n and define its spectral 
distribution function by 

1 n 

J*«( x ) = - Vl(Aj <x), 
n 

i=l 

where t(B) denotes the indicator of an event B. We set F x,I (x) := EJ rXn (x). 
Let g{x) and G(x) denote the density and the distribution function of the 
standard semicircle law 

i rX 

g (x) = — \/4 - x 2 l(\x\ < 2), G(x)= g(u)du. 

For matrices with independent identically distributed (i.i.d.) elements, which 
have moments of all orders, Wigner proved in |10| that F n converges to G(x), 
later on called "Wigner 's semicircle law". The result has been extended in 
various aspects, i.e. by Arnold in [2]. In the non-i.i.d. case Pastur, [9j, showed 
that Lindeberg's condition is sufficient for the convergence. In (7j Gotze and 
Tikhomirov proved the semicircle law for matrices satisfying martingale-type 
conditions for the entries. 

In the majority of previous papers it has been assumed that afj are equal 
for all 1 < i < j < n. Recently Erdos, Yau and Yin and al. study ensem- 
bles of symmetric random matrices with independent elements which satisfy 
n_1 YTj=\ a ij = 1 for all 1 < i < n. See for example the survey of results in jSj. 

In this paper we study a class of random matrices with dependent entries and 
show that limiting distribution for _F Xn (x) is given by Wigner's semicircle law. 
We do not assume that the variances are equal. 

Introduce the cr-algebras 

^ := a{X M :l<k<l<n,(k,l)^ l<i<j<n. 
For any r > we introduce Lindeberg's ratio for random matrices as 

1 n 

L n {r) := ^ ^2E\Xij\ 2 l(\Xij\ > ry/n). 



SEMICIRCLE LAW 



We assume that the following conditions hold 
(1.1) E{Xij\^) = 0; 



(1.2) 
(1.3) 



1 n 

-2 E I ^(Xfj\$ {l,J >) - -> as n -> oo; 
for any fixed r > L n (r) — > as n — )• oo. 



Furthermore, we will use condition (1.3) not only for the matrix X n , but for 
other matrices as well, replacing Xij in the definition of Lindeberg's ratio by 
corresponding elements. 



For all 1 < i < n let Bf := \ YTj=\ °ly We need to 

impose additional conditions 



on the variances afj given by 



(1.4) 
(1.5) 



1 n 

-Y\b. 

i=i 

max Bi < C, 

Ki<n 



1| — > as n — > oo; 



where C is some absolute constant. 



Remark. It is easy to see that the conditions (1.4) and (1.5) follow from the 
following condition 

(1.6) max \Bf — l| — > as n 

Ki<n 



oo. 



The main result of the paper is the following theorem 



Theorem 1.1. Let~K n satisfy conditions ( 1.1 ) — ( 1.5 ) . Then 



sup |F x "(x) - G(x)\ ->■ as n -)■ oo. 



Let us fix i,j. It is easy to see that for all (k, I) ^ 

EXijX kl = EEiXijXultfM)) =EX M E{X ij \^) = 0. 

Hence the elements of the matrix X n are uncorrelated. If we additionally as- 
sume that the elements of the matrix X„ are independent random variables 



then conditions (1.1) and (1.2) are automatically satisfied. The following The- 



orem |1.1| follows immediately in the case when the matrix X n has independent 
entries. 



Theorem 1.2. Assume that the elements X{j of the matrix'K n are independent 
for all 1 < i < j < n and EX tj = 0, EX^ 



conditions (1.3)-(1.5). Then 

sup|F x "(x) -G(x)\ 



afj. Assume that X n satisfies 



as n 



oo. 



The following example illustrates that without condition (1.4) convergence to 
Wigner's semicircle law doesn't hold. 
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Figure 1. Spectrum of matrix X n . 



Example. Let X n denote a block matrix 



X r 



A 



B 
D 



where A is m x m symmetric random matrix with Gaussian elements with 
zero mean and unit variance, B is m x (n — m) random matrix with i.i.d. 
Gaussian elements with zero mean and unit variance. Furthermore, let D be 
a (n — m) x (n — m) diagonal matrix with Gaussian random variables on the 
diagonal with zero mean and unit variance. If we set m := n/2 then it is not 



difficult to check that condition (1.4) doesn't hold. We simulated the spectrum 



of the matrix X n and illustrated a limiting distribution in Figure [TJ 



Remark. We conjecture that Theorem (Theorem 1.2 respectively) holds 
without assumption (1.5). 



Define the Levy distance between the distribution functions F\ and F 2 by 

L(F U F 2 ) = inf{e > : F t (x - e) - e < F 2 {x) < F x {x + e) + e}. 

The following theorem formulates the Lindeberg's universality scheme for ran- 
dom matrices. 

Theorem 1.3. Let X n ,Y n denote independent symmetric random matrices 



with E X 



and EXf, 



afj . Suppose that the matrix X n 



satisfies conditions ( 1.1 ) — ( 1.3 ), and the matrix Y n has independent Gaussian 



elements. Additionally assume that for the matrix Y n condition (1.3) holds. 
Then 

L(F x "(x),F Y "(x)) 



as n — > 00. 



In view of Theorem 1.3 to prove Theorem 1.1 it remains to show convergence 
to semicircle law in the Gaussian case. 
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Theorem 1.4. Assume that the entries Yij of the matrix Y n are independent 

■ihiitinn unih F YL- = (1 F. \ 



for alll < i < j < n and have Gaussian distribution with EYy = 0, E Y^- = of--. 



Assume that conditions dl.3[| — ( 1.5 ) are satisfied. Then 



sup |F Y " (x) — G(x)| — ^ as n -> oo. 



For related ensembles of random covariance matrices it is well known that spec- 
tral distribution function of eigenvalues converges to the Marchenko-Pastur law. 
In this case Gotze and Tikhomirov in [6] received similar results to [ 7j| . Recently 
Adamczak, [I], proved the Marchenko-Pastur law for matrices with martingale 
structure. He assumed that the matrix elements have moments of all orders and 



imposed conditions similar to (1.4). Another class of random matrices with de- 
pendent entries was considered in [8] by O'Rourke. In a forthcoming paper we 
prove analogs of these theorems for random covariance matrices. 



L3 
1.4 



The paper organized as follows. In Section 2 we give a proof of Theorem 
using the method of Stieltjes transforms. In Section 3 we prove Theorem 
by the classical moment method. 

Throughout this paper we assume that all random variables are defined on 
common probability space (ft, J 7 , P). Let Tr(A) denote the trace of a matrix A. 
For a vector x = (x±, ■■■,x n ) let ||sc||2 := (2it=i x i) 1 ^ 2 - We denote the operator 
norm of the matrix A by ||A|| := sup ||Ax||2- We will write a < m b if there 

IMIa=l 

is an absolute constant C depending on m only such that a < Cb. 

From now on we shall omit the index n in the notation for random matrices. 



2. Proof of Theorem 11.31 

We denote the Stieltjes transforms of F x and F Y by S x (z) and S Y (z) re- 
spectively. Due to the relations between distribution functions and Stieltjes 
transforms, the statement of Theorem 11.31 will follow from 



(2.1) \S*(z) - S Y (z)\ -> as oo. 
Set 

(2.2) R x (z) := - zI ) and rY (^) : = (^ Y ~ zI 
By definition 

5 x (z) = — TrER x (z) and S Y (z) = — Tr E R Y (z). 

n n 



We divide the proof of (2.1) into the two subsections 



2.1 



and 



2.2 
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Note that we can substitute r in (1.3) by a decreasing sequence r n tending to 
zero such that 



(2.3) 

and linin-s.oo r n \/n = oo. 



L n (r n ) — > as n — > oo. 



2.1. Truncation of random variables. In this section we truncate the el- 
ements of the matrices X and Y. Let us omit the indices X and Y in the 
notations of the resolvent and the Stieltjes transforms. 

Consider some symmetric n x n matrix D. Put X = X + D. Let 

-l 



R 



1 



X-zI 



71 



Lemma 2.1. 



Tr R - Tr R| < -=■ (Tr D 2 



Proof. By the resolvent equation 
(2.4) R = R 



RDR. 



n 



For resolvent matrices we have, for z = u + iv , v > 0, 



(2.5) 



max{||R||,||R||} < -. 



Using (2.4) and (2.5) it is easy to show that 

1 



Tr R - Tr R| 



1 



n 



Tr RDR| < ^(TrD^)2. 



□ 



We split the matrix entries as X = X + X, where X := X 1(|X| < T n y/n) and 
X :=Xt{\X\ > Tn-s/n). Define the matrix X = {Xij}f j=l . Let 



n(z) 



i , \- x « i 

X-zI and S(z) = -ETrR(z) 
n I n 



By Lemma 2.1 



1/2 



\S(z)-S(z)\< 



1 / 1 



^EX?.l(|Xy|>r nV ^) 



u 2 L^(rn). 



From (2.3) we conclude that 

|5(z) — S(z)\ — > as n — > oo. 
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Introduce the normalized random variables Xij = Xy — E(Xy|5 r ^'^) and the 
matrix X = {Xy}" - =1 . Let 



R(z) :-- 



1 _ V 1 _ 1 

X-zI and S(z) = -ETrH(z) 
n I n 



Again by Lemma 2.1 



1/2 



\S(z)-S(z)\< 



1 / 1 



^EX?.l(|Xy|>r n ^) 



In view of (2.3) the right hand side tends to zero as n — > oo. 



Now we show that (1.2) will hold if we replace X by X. For all 1 < i < j < n 

E(Xj.|^'))-4 

= E(X'.|S (i ' j) ) - E(±?,|3^)) - E(X?-|ffW)) + E(X?.|^')) - 4, 



By the triangle inequality and ( |1.2[ ), (2.3) 

-2 E E i E (4i^' } )-4i 



1 n 

< -2 £ E I E(X||^)) - 4| + 2L n (r n ) 



as n — > 00. 



In what follows assume from now on that IXul < T n y/n and \Yij\ < T n \fn. 



2.2. Universality of the spectrum of eigenvalues. To prove (2.1) we will 
use a method introduced in Ji]. Define the matrix Z := Z(i^) := X cos if + 
Y sin ip. It is easy to see that Z(0) = X and Z(vr/2) = Y. Set W := W(ip) := 
n _1 / 2 Z and 

(2.6) K(z,tp) := (W - zl)- 1 . 

Introduce the Stieltjes transform 



1 n 

S(z,<p) := - E[R(z, ip)]i 



n 

i=l 



Note that S(z,0) and S(z,tt/2) are the Stieltjes transforms 5* x (2;) and S Y (z) 
respectively. 

Obviously we have 

(2.7) s{ ^ ) -s^ ) =f;^ iv . 
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To simplify the arguments we will omit arguments in the notations of matrices 
and Stieltjes transforms. We have 



^ n n 

EE 



dW 1 T^^dZa , 

-e,-e. 



dip \/n ' ^-^ dip 3 ' 

l= X j=i 

where we denote by ej the column vector with 1 in position i and zeros in the 
other positions. We may rewrite the integrand in (2.7) in the following way 

(2.8) ^ = _W^R 

dip n dip 

1 EE^S^r* 



n 3/2 Z-^ Z-^ Qin 
i=l j=l 

i 

-Ui ■ 



i=l j=l 



where Ujj = — [R 2 ]^,. 



£y := Zy, := = - Bin tpXij + cos ipYij, 



For all 1 < i < j < n introduce the random variables 

dZjj 
dip 

and the sets of random variables 

:=U«:l<fc<i<n, (k,l) ^ 
Using Taylor's formula one may write 



where has a uniform distribution on [0,1] and is independent of 
Multiplying both sides of the last equation by 4 and taking mathematical 
expectation on both sides we have 

(2.9) E4- % -(£y,£^)) = E^n^O,^) +E4-^^(0,e W) ) 

+E0(i-^)44.^(^ J ^)). 

By independence of Yij and £W we get 

(2.10) EY ijUij {Q,^ ij) ) =EY lj Eu ij (0,C {lj) ) = 0. 

By the properties of conditional expectation and condition ( |1.1[ ) 

(2.11) £1^(0,^) =E^(0,^'))E(Zy|^) = 0. 
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By (2.9), (2.10) and (2.11) we can rewrite (2.8) in the following way 



1 



n- 



= Ai + A 2 . 
It is easy to see that 

A 1 Q Q O 1 Q 

= _ 2 sin2 ^ij + cos ^XijYij - sin ipXijYij + - sin 2ply. 

The random variables are independent of and ^\ Using this fact we 
conclude that 



dui 



dm 



E^^(o,e fo ' ) ) = 4 E ^(°'^' ) )- 



0. 



(2.12) 
(2.13) 

By the properties of conditional mathematical expectation we get 



(2.14) 



E^-^(0,e (y) )=E^(0,^)E(^|»(«)) 



A direct calculation shows that the derivative of Uij = — [R 2 ]jj is equal to 



du, 



'j 



n2 dZ , ' 
R — — R 



+ 



^eteTR]^ + -=^L 2 ej^B]ji + [Re^R 2 ]^ + -= 

n" J Jn Jn J Jn 



[Re.-efR 2 ]^ 



= [R 2 ] u [R] it H — 7= [R 2 ] /j [R]ii H — 7= [R] ft [R 2 ] jz H — 7= [R] 77 [R 2 ] 

n Jn Jn vn 



Using the obvious bound for the spectral norm of the matrix resolvent ||R|| < 
v^ 1 we get 



(2.15) 



dm 



< 



c 



Jnv s 



From (2.12)-(2.14), (2.15) and assumption (1.2) we deduce 



as n — > oo. 



C J\ 

i A ii<^E E i E (^ (i,i) )-4i^o 

It remains to estimate A 2 . We calculate the second derivative of u 



d 2 u,; 



d£ 2 - 



n2 9V D 5V D 

R "77T — RtTT — R 



RtTT — R "77T — R 



Jl 



R fl^ R ^^ R ^ 



Ti +T 2 + T 3 . 



ji 
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Let's expand the term Ti 

r)V fIV 1 

= T 11 + T 12 +T 13 +T 14 , 



R^R^R 



j i 



(2.16) Ti = -2 

where we denote 
Tii = 

T13 = — [R 2 ],-j[R]jj[R]M, T14 

n n~ 

Using again the bound ||R|| < v^ 1 we can show that 



2 2 

— [R ]jj[R]ji[R]ji, T12 = — — [R ]ji[R]jj[R]iii 

2 - 21 ^ ^ t 2| R 2 ] ii [R] <i [R] ii . 



max(|Tn|, |Ti 2 |, |T 13 |, |Ti 4 |) < 



C 



nv 



l ' 



From the expansion ( 2.16 ) and the bounds of Tu,i = 1, 2, 3, 4 we conclude that 

C 



Ti < 



TO" 



Repeating the above arguments one can show that 

C 



max T 2 , T 3 ) < 



TO 1 



Finally we have 



< 



c 



nv 



4 ' 



Using the assumption < T n y/n we deduce the bound 

|A 2 | < 



From these bounds, ( |2.7| ) and ( |2.8| ) it immediately follows that 

\S(z^-)-S(zM<^- 
We may turn r n to zero and conclude the statement of Theorem 1.3 



3. Proof of Theorem 11.41 



We prove the theorem using the moment method. It is easy to see that the 
moments of F Y (x) can be rewritten as normalized traces of powers of Y: 

k 



J 



x k dF Y (x) = E - Tr ( — J=Y 



1 



n 



1 



It is sufficient to prove that 



E - Tr f — Y 

n \ Jn 



x k dG{x) + o k (l). 



for k > 1, where 0fe(l) tends to zero as n — > oo for any fixed k. 
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It is well known that the moments of semicircle law are given by the Catalan 
numbers 



Furthermore we shall use the notations and the definitions from j3j. A graph 
is a triple (E, V, F), where E is the set of edges, V is the set of vertices, and F 
is a function, F : E — > V x V. Let i = be a vector taking values in 

{1, n} k . For a vector i we define a T-graph as follows. Draw a horizontal line 
and plot the numbers ii, on it. Consider the distinct numbers as vertices, 
and draw k edges ej from ij to = 1, k, using = i\ by convention. 

Denote the number of distinct ij's by t. Such a graph is called a T(k, t)-graph. 

Two T(k, i)-graphs are said to be isomorphic if they can be converted each other 
by a permutation of (1, n). By this definition, all T-graphs are classified into 
isomorphism classes. We shall call the T(k, t)-graph canonical if it has the 
following properties: 

1) Its vertex set is {1, ....,£}; 

2) Its edge set is {e±, e^}; 

3) There is a function g from {l,...,k} onto {l,...,t} satisfying g(l) = 1 and 
g(i) < max{<7(l), g(i — 1)} + 1 for 1 < i < k; 

4) F(ei) = (g(i),g(i + 1)), for i = 1, k, with the convention g{k + 1) = <?(1) = 
1. 

It is easy to see that each isomorphism class contains one and only one canonical 
T-graph that is associated with a function g, and a general graph in this class 
can be defined by F(ej) = (i g (j), i g (j+i))- It is easy to see that each isomorphism 
class contains n(n — l)...(n — t + 1) r(fc, t)-graphs. 

We shall classify all canonical graphs into three categories. Category 1 consists 
of all canonical T(k, t)-graphs with the property that each edge is coincident 
with exactly one other edge of opposite direction and the graph of noncoincident 
edges forms a tree. It is easy to see if k is odd then there are no graphs in 
category 1. If k is even, i.e. k = 2m, say, we denote a t)-graph by Ti(2m). 
Category 2 consists of all canonical graphs that have at least one edge with 
odd multiplicity. We shall denote this category by T2(k,t). Finally, category 3 
consists of all other canonical graphs, which we denote by T3(k,t). 

It is known, see (3l Lemma 2.4], that the number of Fi(2m)-graphs is equal to 

1 /2m\ 
m+l \ m J ' 

We expand the normalized traces of powers of Y in a sum 



where the summation is taken over all sequences i = i^) S {1, n}. 




(3.1) 
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For each vector i we construct a graph G(i) as above. We denote by Y(i) = 
Y{G(i)). 

Then we may split the moments of -F Y (x) into three terms 

E k Tr (i Y ) fc = ni^^^"^ 1 =S 1 + S 2 + S 3 , 
\ v / ^ 

where 

^■ = ^+1 E E 

r(fc,t)ec j G(i)&(k,t) 

and the summation ^r(fct)ec * s taken over all canonical T(k, i)-graphs in cat- 
egory Cj and the summation ^G(i)er(fe t) ls taken over all isomorphic graphs 
for a given canonical graph. 

From the independence of Y^ and E Y^ s_1 = 0, s > 1, it follows that S% = 0. 

For the graphs from categories C\ and C3 we introduce further notations. Let 
us consider the T(k, t)-graph G(i). Without loss of generality we assume that 
ii,l = l,...,t are distinct coordinates of the vector i and define a vector it = 
(it,..., i t ). We also set G(i t ) := G(i). Let i t = (ii, i q -i, i q+ i, it) and 
it = (ii, ip-i, ip+i, ig-i, ig+i, it) be vectors derived from i t by deleting 
the elements in the position q and p, q respectively. We denote the graph without 
the vertex i q and all edges linked to it by G{%t). If the vertex i q is incident to 
a loop we denote by G'{it) the graph with this loop removed. 

Let it denote a vector taking values in {1, ...,n} 1 and j t and j t denote vectors 
from {1, n}* -1 and {1, n}* -2 respectively, derived from j t by deleting the 
elements in the position q and p, q respectively. By a graph G(j t ) (respectively 
G(j t )) we mean the graph G(it) (respectively G(it)), while allowing repeated 
vertices. 

Now we will estimate the term S3. For a graph from category C3 we know that 
k has to be even, i.e. k = 2m, say. We illustrate the example of a r3(fc, t)-graph 
in Figure[2j This graph corresponds to the term Y(G(iz)) = Y^Y^Y^Y^. 
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We mention that EK 2 ? <« of s - . Hence we may rewrite the terms which 
correspond to the graphs from category C3 via variances. 

In each graph from category C3 there is at least one vertex incident to a loop 
with multiplicity greater or equal to two or an edge with multiplicity greater 
then two. It is possible as well that both cases occur. 

Suppose that there is a vertex, let's say i\, which is incident to a loop with 
multiplicity s > 2. It remains to consider the remaining 2(m — s) edges. We 
will consequently delete edges and vertices from the graph using the following 
algorithm: 



(1) If the number of distinct vertices is equal to one we should go to step 
(3). Otherwise, we take a vertex, let's say i q ,q / 1, such that there 
are no new edges starting from it. We also take a vertex i p which is 
connected to i q by an edge of this graph. There are three possibilities: 
(a) There is a loop incident to vertex i q with multiplicity 2a, a > 1. In 
this case we estimate 

-Lj £ e[y(g(i))] < a ^E^n^Kl- 

G(i)eT(2m,t) it 



Applying the inequality n l af i <Bf <C 2 ,a times we delete all 
loops incident to this vertex; 
(b) There is no loop incident to i q , but the multiplicity of the edge 
from i p to i q is equal to 2. In this case we estimate 



1 n 

E nY(G(i))] <^Ew t ))iE<v 



n 

G(*)er(2m,t) 



We may delete the vertex i q and the two coinciding edges from i p 
to i q using condition (1.5); 
(c) There is no loop incident to i q , but the multiplicity of the edge 
from i p to i q is equal to 26, b > 1. In this case we estimate 



G(i)er(2m,t) i t iq=l 



n m+l n 



Here we may use the inequality n l (jf p i q < Bf^ < C 2 , 6 — 1 times 
and consequently delete all coinciding edges except two. Then we 
may apply (b); 
(2) go to step (1); 
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(3) in this step one may use the bound 

i n i n n2(s-l) n 

1 \ " w v 2s ^ 1 \ ~* 2s <- ° \ " 2 

«i=l ii=l ii=l 

r»2(s-i) n 
n 11=1 

<C 2 (- 1 )r„ 2 + C 2 (- 1 )L n (r n ) = 0s (l), 
where we have used the inequality n^af , < B 2 ^ < C 2 . 

Applying this algorithm we get the bound 

£ E[7(G(i))] < m C 2 ( m - 1 )(r n 2 + L n (r n )) = 0m (l). 
G(i)er(2m,t) 

If there are no loops, but just an edge with multiplicity greater then two, then 
we can apply the above procedure again and use in the step (3) the following 
bound 



I n l n (72(s-2) n 

S+T Y E ^i2 ^ — h Yl ^ „3 Y ^2 



11,12=1 11,12=1 «1,«2 = 1 

il^«2 



81,12=1 11,12=1 

<C 2 (- 1 )r 2 + C 2 ( s - 1 )L n (r n ) = o s (l), 

where we have used the inequality n a 2 , < 5 2 < C 2 and (1.5). 

As an example we recommend to check this algorithm for the graph in Figure [2j 

It is easy to see that the number of different canonical graphs in C3 is of order 
O m (l). Finally for the term £3 we get 

S3 = o m (l). 



It remains to consider the term S\ . For a graph from category C\ we know that 
k has to be even, i.e. k = 2m, say. In the category C\ using the notations of 
it,it,h,jt,3t and j t we take t = m + l. 

We illustrate on the left part of Figure [3] an example of the tree of noncoincident 
edges of a ri(2m)-graph for m = 5. The term corresponding to this tree is 
YfCda)) = Y 2 Y 2 Y 2 Y 2 Y 2 

We denote by a 2 (i m+ i) = a 2 {G{i m+ \)) the product of m numbers aj i , where 
i s ,h,s < t are vertices of the graph G(i m+ i) connected by edges of this graph. 
In our example, a 2 (i m+1 ) = a 2 (i 6 ) = cr 2 ii2 a 2 2i3 a 2 2 i4 ^ 2 i 5 4 V 



SEMICIRCLE LAW 



15 



Figure 3. On the left, the tree of noncoincident edges of a 
ri(10)-graph is shown. On the right, the tree of noncoincident 
edges of a ri(10)-graph with deleted leaf iq is shown . 



If m = 1 then a 2 'fa) = (T 2 , and 

( 3 - 2 ) ^ E a hi2 = - E - E a li2 - 1 



11,12=1 



n 1 — ' n 

ii=l L 12=1 



+ l + o(l), 



where we have used n~ 2 =1 a hii = °(1)- By (1-4) the first term is of order 
o(l). The number of canonical graphs in Ci for m = 1 is equal to 1. We 
conclude for m = 1 that 

n 

5 i = n ^ 2 E E ^2 = l + °( 1 )' 

Ti(2) »i,*2=l 

n^«2 

Now we assume that m > 1. Then we can find a leaf in the tree, let's say 
i q , and a vertex j p , which is connected to i q by an edge of this tree. We have 

cr 2 (i m +i) = cr 2 (im+i) ■ crf piq , where o 2 (i m +i) = a 2 (G(i m +i))- 

In our example we can take the leaf i§. On the right part of Figure [3] we 
have drawn the tree with deleted leaf ia. We have of { = a 2 , ■ and a 2 fa) = 



2 2 2 2 
a i\ii a i<ziz° iiiA U 



It is easy to see that 

( 3 - 3 ) ^+1 E °" 2 (Wl) = "7^1 E ff2 (im+l) + °m(l) 



Jm + l 



7? 



1 - 

^+1 E ^Om+l) E + 



3m + l 



— y o 2 q 



m+li 



1 ™ 

- V °f i - 1 



3m + l 



i„=l 



+ ^E ct2 (^ 



t+l) + Om(l), 



Jm + l 
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where we have added some graphs from category C3. Now we will show that 
the first term in (3.3) is of order o m (1). Note that 



(3.4) 



E 0-2 Cwi) 



3m + l 



n 

-E 



1 H 



n 



1 n 

-Y oh -1 



n 



—[ E °" 2 (jm+l)- 



Jm+1 



We can sequentially delete leafs from the tree and using (1.5) write the bound 
(3-5) ~ £ -^+1) * 



3m+l 



By (3.5) and (1.4) we have shown that (3.4) is of order o m (l). For the second 
term in (3.3), i.e. 



— E 0-2 Owi) 

Jm+l 

we can repeat the above procedure and stop if we arrive at only two vertices in 

Fina 



the tree. It the last step we can use the result (3.2). Finally we get 

^TT E £^ 2 (Wi) 



Si 



F 1 (2m) im+i 



1 



m + 1 V tn 



+ o m (l), 



which proves Theorem 1.4 
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